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Abstract 

We reinforce the observations of almost stable scattering in noninte- 
grable models and show that PT-symmetry can be used as a guiding 
principle to select relevant systems also when it comes to integrability 
properties. We show that the presence of unbroken TT-symmetry in 
classical field theories produces quasi-integrable excitations with asymp¬ 
totically conserved charges. 


1 Introduction 

Since it was first conjectured [T], the use of TT-symmetry as a guiding 
principle to identify potentially relevant physical systems has been applied to 
various quantum settings but there have been also investigations on the effect 
such a symmetry has at the classical level. 

A vast range of classical models was studied, from point particle systems 
[2] to classical field theories [SI IH 13 El [Z] • An interesting observation was that 
the requirement of symmetry under parity and time reversal of Calogero-like 
particles would imply that a class of solitons of the Boussinesq equation would 
have real energies [8]. 

At that time it was noticed that most classical field theories of interest 
were described by equations which were TT-symmetric. There have been many 
relevant works investigating TT-symmetric integrable models, e.g. liiiniiiii 
da [13 [H dsi na [m [18]. However, this property has not been scrutinised 
thoroughly enough and there are aspects which remain little explored. For 
instance the transition between the broken and the broken phases of the VT- 
symmetry has not been fully understood. 

In a quantum system, it was shown [I] that the invariance of the Hamiltonian 
with respect to parity and time reversal is not enough to guarantee the reality 
of energies. In order to have real eigenvalues it is necessary to have invariant 
eigenstates as well. If this is the case, eigenenergies are real, otherwise they 
would come in complex conjugate pairs. 
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Classically, what are the effects of being in the unbroken or broken phases 
of a T^T-symmetric model? It turns out this symmetry can impose restrictions 
on the phase space which can mimic those of integrable systems giving rise to 
almost stable solitonic solutions. Quasi-integrable models have been studied 
independently, as in m, since a particular generalisation of the sine-Gordon 
model was proposed |20] . Initially just a numerical observation, the appearance 
of quasi-integrable behaviour has been more carefully investigated and a deeper 
comprehension of the reasons why charges were almost conserved was pursued. 
This led to an interesting algebraic description of such a phenomenon in [T5] . 
where the authors identified the intriguing role played by VT. 

Here we show how one can understand the quasi-integrability of unbroken 
T^T-symmetric problems in terms of the flat gauge connections in the Lax for¬ 
mulation of intgrable models. In the next section, we motivate the present 
investigation by showcasing a series of well known integrable nonlinear partial 
differential equations with a few important solutions and call attention to the 
space-time symmetries of those. Then, in section [3] we review the construction 
of conserved charges based on a path-independent solution of Lax’s auxiliary 
problem. These tools are used in section |4l where a particular nonintegrable 
model is investigated and shown to possess almost PT-symmetric excitations 
which lead to the quasi-conservation of integrals of motion. In section [S] we 
construct the approximate conserved charges based on the asymptotic validity 
of the inverse scattering method and finally, in section [51 a numerical procedure 
is employed to check the approximate constancy of the spectral determinant. 
We finish by stating our final remarks in section [T] 

2 PT-symmetry in integrable field theories 

As we have argued, it is know for years that important models of classical 
field theory are associated to invariant equations under time reversal and parity 
operations. In what follows we revisit some of them and highlight their be¬ 
haviour under VT. More than that, we show that their most relevant solutions, 
namely their solitonic excitations, are also invariant. Thus, the regime of stable 
particle-like waves which can scatter indefinitely, preserving the form of their 
profiles and conserving energy, coincide with the so called unbroken phase of 
VT -symmetry. 

As introduced in [1], the T^T-symmetry of a Hamiltonian system. 


[H,vr] = o, 


( 1 ) 


is not enough to guarantee the reality of its spectrum; instead one should also 
impose some restriction on the eigenvectors. If the eigenstates are also invariant, 
up to a phase, under this symmetry. 


( 2 ) 


VTip = ±ip 


one is said to be in the unbroken PT-symmetry regime. It turns out that 
this concept can be transported to ensure that the energy of a specific field 
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configuration in classical mechanics is real [J]. The analogous expression for the 
unbroken VT symmetry condition is then 

VTu = ±u. (3) 

Our findings suggest that some systems that are not exactly integrable but 
which satisfy an invariance of this type can behave as quasi-integrable, in the 
sense that they possess families of asympotically preserved charges. Before 
studying the behaviour of PT-symmetry with respect to deviations from the 
integrable domain, we present in the next subsections completely integrable 
systems and highlight their space-time invariances. 


2.1 Korteweg - de Vries 

The KdV equation has been well studied also in the realm of non-Hermitian 
Hamiltonians. In [5], it was shown that a certain PT-symmetric deformation 
admits the existence of certain conservation laws and in [4] it has emerged that 
a family of Hamiltonian flows can be associated to a particular deformation of 
the KdV equation. Further generalisations were considered in and in [7] 
compactons were investigated. 

But one does not have to impose PT-symmetric extensions of the KdV in 
order to appreciate the effects of these simultaneous space-time transformations. 
In fact, taking the standard equation. 


Ut CLUUxx ^^xxx — Hi 

arising form the Hamiltonian density 


a 

—1 

6 


TL = --u^ + -ui, 


( 4 ) 

( 5 ) 


with usually o = — 6 and 6 = 1, one can see that the Hamiltonian remains 
ivariant if we replace x —> —a;, t—> —t and u —> m at once. 

This model admits multi-solitons excitations, and a two-soliton solution is 
known to be 


m(x, t) 


_12^^^ (c 2 sinh^ (^i) -b ci cosh^ fe))_ 

ii'/f - \/f-) cosh (Cl -b 6 ) + {^/f + ) cosh (Cl - C2)) 2 ’ 


where we have abbreviated Ci = ~ Cit). The action of the symmetry 

operation into these objects gives VT[^i] = —Ci making this solitonic solution 
u{x,t), which is even in Ci, invariant, i.e. VT[u{x,t)] = u(—x, —t) = u{x,t), as 
required. In Figure [Done can see the space-time symmetry of the two-solitons 
solutions. 

These properties can be immediately identified if one plots the evolution 
surface for this KdV solitons collision, which we show bellow. There is evidently 
an axis of symmetry with respect to which the simultaneous reflections on the 
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Figure 1: Colision of KdV solitons, repulsive (bounce-exchange) on the left 
and attractive (merge-split) on the right. The solution is PT-symmetric with 
respect to the centre of mass reference frame. 


space and time coordinates produce an identity operation, be it a repulsion or 
an attraction. 

Like the KdV model, there are many other examples of integrable field equa¬ 
tions which not only present the parity-time reversal symmetry bu also admit 
T^T-symmetric stable solutions. In what follows we discuss a few more illustra¬ 
tions. 


2.2 Sinh-Gordon 


The Sinh-Gordon models is described by a Lorentz-invariant equation, 


4>tt - 


+ ^ sinh/?(/) = 0, 


(7) 


derived from 

7^ = ^ (</>?+ </'x) + ^ (cosh/30 - 1). 


( 8 ) 


In the latter, the derivative terms are invariant under space-time reflections and 
all terms remain unchanged if solutions belong to the unbroken PT-symmetry 
phase, PT0 = ±0. Remarkable solutions for this equation present a travelling 
wave form, stable under collision, as depicted bellow. Some of these complex 
waves can be cast in the following form 


(j){x,t,UJ,fx) = -log 


I — 2 COs(w)e^^* _j_ g2/i(tsinh(e)-|-£!;cosh(e)) \ 
^ _j_ 2 (;;Qg(^)gA^(i sinh(^)+a: cosh(^)) _j_ g2p(i sinh{0)-t-a: cosh(^)) J 


+ 


2i 

— arctan 


2 

g2/j(i sinh(0)-t-a: cosh(0)) _ 


(9) 


so the PT-invariance of the equation of motion above, as a: —>■ —x and t —>■ —t, 
reads 


(j){x,t,ui,fj,) = <j){-x,-t,-u},fj,) = <j){-x,-t,uj,fi)*, (10) 
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in accordance with the complex conjugation that comes along with time-reversal 
in quantum, and complex, field theories (see Figure[2]for the real and imaginary 
parts of the excitation). Here we have used 9 = tanh“^ (^) as the rapidity 
variable, with c = 1 . 




Figure 2: Sinh-Gordon complex solution (real part on the left, imaginary part 
on the right). 


2.3 Sine-Gordon 


By complexifying the coupling constant /3 —>■ *7 of the previous model, one 
obtains the Sine-Gordon Hamiltonian density, 

1 777 ^ 


and its equation of motion, 


(t>tt 


no . , ^ 

-\ -sinycp = 1 ), 
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with localized solutions, 




( 12 ) 


(13) 


describing kinks (e = -1-1) and anti-kinks (e = —1), depicted in Figure [31 In 
order to appreciate its VT properties one should note that 


4 (arctan (e^) — arctan (—e '")) = 0 mod 27r. 


(14) 


An interesting feature is therefore that these two fundamental configura¬ 
tions, the kink and anti-kink, from which one can construct other excitations 
(breathers, bound states), are connected via T^T, up to a topological equiva¬ 
lence. The system composed of a single (anti)-kink is not explicitly PT-invariant 
but if one wants to investigate their stability one can create a collision between 
a kink and an anti-kink and this process is similar to the one shown in Figure 
[2] (on the right). Other non-symmetric collusions are also possible but in those 
cases topological arguments are to be used. 
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Figure 3: Sine-Gordon kink (dashed line, in black) and anti-kink (continuous 
line, in gray), interpolating neighbouring vacua. 


2.4 Boussinesq equation and Calogero scattering 

Finally, we illustrate the good VT properties of integrable models by pre¬ 
senting one last example, that of the Boussinesq equation. In [8] we have ex¬ 
plicitly constructed a class of solutions depending on dynamical movable poles. 
Interestingly, the poles satisfy an interaction dictated by the Calogero Hamil¬ 
tonian m- In that occasion we have concluded that PT-symmetric Calogero 
particles would guarantee that the solutions to the Boussinesq equation have 
real energies. 

Below we depict the evolution of the real and imaginary parts of the Boussi¬ 
nesq PT-symmetric solution associated to the existence of three Calogero-like 
poles. 


JV/\— 
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Figure 4: Rational Boussinesq solitonic solution with Calogero poles: real (left) 
and imaginary (right) components, [8]. 


One can verify from Figure 0] that the fields are indeed invariant under 
space and time reversal: the real part at two points t and —t correspond to a 
simple reflexion x —>■ —x, but the imaginary part also suffers a — sign inversion, 
associated to the complex conjugation. 

The poles appearing in this solution interact as shown in Figure [5l If we 
decompose the scattering in two-bodies collisions, the trajectories of the particles 
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are symmetric with respect to space-time inversions relative to the centre of mass 
reference frame. 



Figure 5: Trajectories of a Calogero three-particle system in a classical scatter¬ 
ing, repulsive (bounce-exchange) on the left and attractive (merge-split) on the 
right, [H]. 


Although constructed from a classical starting point, these ideas can be ex¬ 
tended to quantum scattering processes. The classical collisions above resemble 
the decomposition property of the scattering matrix. From Figure [S] we see the 
effect of space-time reversal and it becomes evident that indeed VT^ = 1. 



Figure 6: Two-particle scattering invariant under V'T transformation, up to a 
minor relabeling. 


In the case of quasi-integrable systems, despite the asymptotic stability of 
the excitations one observes the presence of radiation near the interaction re¬ 
gion. It has been noticed inSEillEl that unstable particles can be produced 
and annihilated without being observed asymptotically in a scattering process. 
These particles are associated to complex poles in the scattering amplitude and 
could be related to the present observations. 
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3 Path-independence and conservation laws 

Once appreciated the role of PT-symmetry in integrable systems, we now 
turn our attention to an attempt to understand the origin of important associ¬ 
ated properties, more notably the stability of solutions. In order to do so, we 
must start from a formulation of integrable systems in which the space and time 
coordinates arise naturally. In fact, an integrable theory can be formulated in 
terms of Wilson loops in space-time and the non abelian Stokes theorem , 


W{T) = Pexp 

— ® Andx^ 

= V exp 

- [ dx’^dx>^W-^Ff,^W 


Jr 


Jt. 


= ns). (15) 


One can have closed paths starting (and finishing) at different points. For 
instance, using the four vertices - (—F, — r), {+L, —r), {+L, -l-r), (—L, -|-t) - of 
Figure [71 we can write, for a closed contour around the point (—L, —t). 


W(Fc_) = W{T_L)W{T+r)W{T+L)W{T_r). (16) 


t 


+ T 


W{T+r) 


-L 


iF(r_i) 




W{T^r) 


PT-symmetry axis 



W{V+l) 


-\-L 


X 


Figure 7: A closed path will in general produce a phase shift in (1251) due to the 
presence of the potentials in that region of space. However, if these gauge 
potentials present certain behaviour under parity reflection and time reversal 
then the integral in (1261) can be path-independent. 












We note that the object W in (fT3|) will depend on 


J dtAt{L,t) + j dtAt{—L,t) = j dtAt{L,t) — J dtAt{—L,—t) (17) 


and 


J dxAx{x,-T)+ J dxAx{x,T) = — J dxAj;{-x,-T) + J dxAx{x,T) (18) 


which will vanish if these gauge potentials are invariant under parity and time 
reversal operations, 


^t{L, t) — At{—L, —t). 

(19) 

Axix,T) = At{-X,-T). 

(20) 

Indeed, then we have, respectively, 


U{t) = W(T+l) = W(T_l)-\ 

(21) 

Six) = W(r+,) = wiT_r)-\ 

(22) 

and consequently 


WiTc_) = Uit)-^Six)Uit)Six)-\ 

(23) 

can be associted to a conserved object, since 


detW(rc_) = 1. 

(24) 

If we have an auxiliary vector ?/; of the following form 


-ipixA) = fF(r)^/>o, 

(25) 


with 


W(r) = 'Pexp^y dxdtW{x,t) ^Fxt{xA)W{xA^ , (26) 

where the nonintegrable phase of this path-ordered exponential is dictated by 
the object Fxt{xA)^ then its value will heavily rely on the chosen curve F. If 
Fxt{x,t) vanishes, the system becomes integrable as the phase trivialises and 
the auxiliary vector satisfies 

[dt +At{x,t)]'il^{x,t) = 0, (27) 

[dx+Ax{x,t)]il^(.x,t) = 0, (28) 

and the object Fxt(x,t) can be identified with 

Fxt = [dt + At^dx + Ax\ = 0 . 
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(29) 


If one is concerned with closed paths F c there is another way of obtaining a 
trivial integrated phase, W{Tc) = 1) and that is if the integrand is symmetric 
with respect to space and time along a certain integration contour, mimicking 
the effects of a zero curvature condition F^t = 0. That means that, albeit not 
being necessarily integrable, a model might present excitation configurations 
that resemble in some senses those of an integrable theory. 

So far we have presented various examples of integrable systems, real and 
complex, and we have highlighted their PT-symmetric properties. In what 
follows we somehow reverse the logic and present some examples of nonintegrable 
models and argue that if PT-symmetry is present then we can appreciate the 
survival of a few features which are reminiscent of integrable systems. 

To achieve this we have chosen a deformation of the sine-Gordon model 
which was introduced in |20) and further explored in |19| . where stability and 
symmetry properties have been spotted. 


4 The PT-integrabitliy of a deformed sine-Gordon 

Taking specifically the example of the sine-Gordon model, we can write the 
gauge fields associated to the laboratory coordinates as in [TH] . For our purposes 
here, in investigating the Inverse Scattering problem for quasi-integrable models, 
it will be more convenient to express the problem in terms of space and time 
coordinates, rather than light-cone coordinates. Starting from a standard form 
of a more symmetric Lax pair associated with the sine-Gordon model, we express 
it in terms of a generic potential The result are the following gauge fields, 

+ ia+, (30) 

2i 

Ax{x,t) = +ia-, (31) 


with 


a± 


1 

2y/2 



v'mx,t)) _ _ ^ 
\/V{^{x,t)) 


Vvmx,t))c 


,(32) 


making the potentials invariant under the operations of space-time reversal, 
m and (EOI), if so are the fields. 

As a matter of fact, these properties also render the equations (HZl) and (1281) 
T^T-symmetric if one considers the parity transformation to act on the 2x2 

Pauli matrices as V = 



VaAP = <Jx, TaAT 

VcTyV = —CTy, TcTyT 

Va^V = -cr^, TctzT 


t 

~(Xy, 

(Xz, 


VTa^VT = dx, 

VTdyVT = CTy, (33) 

VTdAPT = 


10 







Figure 8: Solitary-like waves collision for n=2,3,4,5 depicted in terms of the 
fields $ (left) and momenta If (right). As the deformation parameter e = n —2 
increases, the less stable the scattering tends to be and some radiation starts to 
be observed after the collision, breaking the time reversal symmetry. There is 
also a slightly better behaviour of positive n compared to odd values. 


The unbroken VT -symmetry of the classical field is thus not directly related to 
the unbroken phase of the vector tp, but with the PT-symmetry of the equation 
itself. 

The form of the potentials in (l30|) - (IM]l is convenient as it allows one to 


11 

























































































































treat nonintegrable extensions, such as [20] 


2 




32 (4M)2 

^ (2S)2 


tan^ 




(34) 


depending on a deformation parameter e = n—2, which vanishes for the the sine- 
Gordon case, with n = 2. Such gauge connections correspond to an anomalous 
curvature, 


Fa:t = + y'($(a;,t))) cr^ + A, 


(35) 


A = 


X 


2V2B V 2 V 


2V''{^{x,t)) l/'($(a;,t))2 \ 

vvmx,t)) ~ vmx,t))3/^) 


(36) 


It is interesting to note that, again, the zero curvature condition depends 
not only on the model - as it is the case for integrable systems - but also on the 
solution. However, in order to have a vanishing anomaly one needs constant 
solutions, $t(a;,t) = ^x{x,t) = 0, which necessarily occurs asymptotically for 
models with degenerate vacua. Unless, of course, if U = Vsg, when the Lax 
pairs are those commonly found in the literature, e.g. [53], and the second factor 
of A is identically satisfied. 

From Figure [5] indicate that as the deformation parameter increases there 
is a tendency of forming radiation waves, as observed in |19j . These excita¬ 
tions cannot be identified with unstable particles nil Hi HZ], as the latter are 
present in integrable systems and the former constitute a sign of the breaking 
of integrability. 


5 Symmetries of spectral properties 

An important property of the spectral coefficients a(A) and 6(A) is that 
the action-angle variables are constructed in terms of them. This implies one 
can expand the action variable in powers of the spectral parameter A, so that 
if action variable is constant then one can generate infinitely-many conserved 
charges. Relaxing this latter condition a bit, if the action variable changes 
slowly with time keeping its final state equal to its initial state, then we have 
an infinite number of quantities whose values are asymptotically the same. 

Thus, if one is concerned only with the in- and out-states in a far past and 
a far future, respectively, there is the possibility of constructing a great deal of 
objects which are conserved in a broader sense. In order to so, it is essential to 
investigate the behaviour of the spectral data, a(A) and 6(A). We can expand 
at a; —>■ 00 


'ip+{x,A) = a{A)tj;-{x,A) + b{A)ip-{x, A), (37) 
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Figure 9: Solitary waves scattering for the potential in equation (1341) with n = 
2,3,4, 5 clockwise from top left. In the vertical axes we have plot the space 
coordinates (x) and in the horizontal line the time coordinates (t) flow; the 
amplitude of the wave at each point is representated in greyscale, for which the 
larger amplitudes are plotted in darker shades. 



Figure 10: Diagonal components of Ax{x, t) for interaction potential in equation 
(|M1) with n = 2, 3,4, 5 clockwise from top left. Similar figures are produced 
when computing the off-diagonal components of Ax{x,t), as well as diagonal 
and off-diagonal components of At{x,t) 


and apply to it the time evolution. 

As a consequence, one can show that, if the phases tend to multiple values 
of the asymptotic vacua topological charge, then the left hand side of the first 
equation identically vanishes and we automatically have that the a (A) keeps its 
original value, ao(A), if one is interested in asymptotic states only. On the other 
hand, as can be seen in more detail in [23] , the other spectral coefficient evolves 
according to the potential I4(<I>), evaluated at the initial vacua configurations 


a(A, t) = ao(A), &(A, t) = 60 (A) (38) 
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with an oscillating phase, 

uj(A) r($) 


(39) 


4V2M 

and an exponential decay/growth term, 



(40) 


The last terms do not contribute if the minimum of the potential, achieved 
by the field asymptotically, can be set to vanish by a simple shift - it is worth 
reminding that, even in this situation, the first term remains finite. Then the 
r factor does not affect the time evolution and we only have the oscillating 
behaviour left. Since the conserved charges can be obtained as the coefficients 
in a power series of the spectral parameter A for a function of a(A) only, then one 
can argue that the charges themselves must also be asymptotically conserved in 
time. 

In fact, a simple numerical evaluation shows that, for the model proposed 
by Bazeia et al, the spectral curve at different instants indicate a small non¬ 
trivial time evolution of the scattering data. However, when we compare the 
asymptotic states the spectral curves tend to coincide, as seen in Figure ITT] 

6 Numerical validation 

In order to verify the results previously stated, here we present some nu¬ 
merical analysis to investigate the scattering of solitary waves which present 
approximately PT-symmetry. In Figure |9] we show a density plot of the field 
excitations as time evolves along the horizontal axis, as colliding wave pack¬ 
ets seen from above. The dynamics is governed by a generalised sine-Gordon 
interaction of the type found in [20] . 

Having the excitations an approximate P7~-symmetry we can check for the 
T^T-invariance on the gauge potentials Ax(x, t) and At{x, t). Therefore, we show 
in Figures uni ??, ??, ?? the behaviour of the diagonal and off-diagonal com¬ 
ponents of theses connections and confirm the desired properties are observed, 
at least in a first approximation. 

In fact, it becomes evident, however, that as the deformation parameter, 
e = n — 2, moves away form the integrable situation, n = 2, the scattering 
symmetry tends to deteriorate. 

In all of these cases, one could also plot the behaviour of the anomaly A 
coming from the nonflatness of the curvature and it is possible to verify 
that its values is relatively small and bounded in this region we are interested. 

Finally we present the spectral curve for the model proposed by Bazeia et 
al [20], at two different instants. The figures show that there is indeed a smalll 
deviation from the conserving regime but the spectral curve, and consequently 
the conserved charges, tend to recover their assimptotic values. 
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Figure 11: On top we show the evolution of an initial configuration under the 
Bazeia potential for n = 4,6, respectively. Then, we depict the spectral curves 
at five different instants, equally separated, represented by different colours. 
We can observe it is not constant in time but one can see after some time the 
spectral curve tends to recover the initial profile. Finally, at the bottom the 
spectral curves at approximately symmetric asymptoitic times. The greater the 
deviation from the integrable case (n=2) the more the spectral curves vary. 


7 Concluding remarks 

The space and time coordinates play a key role in the zero curvature for¬ 
mulation based on the nonabelian Stokes theorem. In this work it was shown 
that investigating the VT symmetry of an integrable model at that level might 
prove interesting as one considers deformations which in principle break the in- 
tegrability of the system. We have presented a general scheme and a particular 
example where asymptotic charges may be considered to be preserved even if 
the model does not admit a compatible Lax pair. The numerical studies car¬ 
ried out also emphasize the connection between PT-symmetry and the so called 
quasi-integrahility. 
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